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Exact droplet-like solutions to the nonlinear scalar field equations have been obtained 
in the Robert son- Walker space-time and their linearized stability has been proved. 



1 Introduction 



When the general theory of relativity (GTR) and quantum theory of field were devel- 
oped, an interest to study the role of gravitational interaction in elementary particle 
physics arose. On this context, to obtain and study the particle- like solutions to the 
consistent systems of wave and gravitational fields present a major interest. To obtain 
and study the properties of regular localized solutions to the nonlinear classical field 
equations (soliton- or particle-like solutions) is connected with the hope to develop a 
divergence-free theory of elementary particle, which in its turn would describe the com- 
plex spatial structure of particle, observed experimentally. In doing so one should keep 
in mind that the nonlinear generalization of field theory is necessary irrespective of the 
question of divergence as the consideration of interaction between the fields inevitably 
leads to the advent of nonlinear terms in the field equations. Consequently, nonlinearity 
should be considered not only as one of the ways to eliminate difficulties of theory, but 
also the reflection of objective properties of field. As it is noticed by N. N. Bogoluibov 
and D. B. Shirkov [1], the complete description of elementary particles with all their 
physical characteristics (say, magnetic momentum) can give only the interacting field 
theory. So one can say that individual free (linear) fields present themselves as the basis 
to describe these particles in the framework of interacting field theory. As elementary 
particle is a quantum object, so the attempts to develop a classical model of particle 
remain preliminary but necessary stage of study for transformation to quantum theory. 

In this paper the system of interacting scalar and electromagnetic fields are being 
considered in the Robertson- Walker Universe with the metric [2] : 



ds 2 = dt 2 - R 2 {t) 



dr 



2 



1 — kr 2 



+ r 2 d6 2 + r 2 sin 2 #d0 2 



(1.1) 



where R(t) defines the size of the Universe, and k takes the values and ±1. Droplet: 
it is some kind of soliton-like solutions to the field equations possessing sharp boundary. 
Similar solution was first obtained by Werle [3]. Further, a series of work was done 
where the solutions with sharp boundary to the nonlinear field equations were being 
found and studied in external gravitational field as well as in the selfconsistent one 
[4-10]. Present paper generalizes the partial results obtained by the authors earlier. 
Moreover here the question of stability is considered which presents a growing interest. 

2 Fundamental equations and their solutions 

We will choose the Lagrangian of interacting scalar (</?) and electromagnetic {J- a p) 
fields in the form [4]: 

C = (l/2)wv?' Q - (1/4)^ (2.1) 

where the function ^(tp) = 1 + K<S>(tp) characterizes the interaction (^(y?) = 1 cor- 
responds to the system of free fields). We will seek the static spherically symmetric 
solutions assuming that the scalar field tp is the function of r only, and the vector field 
Afx possesses one component Ao(r), i.e. 

if = tptr), A, = 6°A (r) = %A(r). 



2 



It means that (Tap) also possesses one component i.e. 

= (5l5l - 5l5l)T 01 (r) = A!(r), 

where / denotes differentiation with respect to r. 

The equations to scalar and electromagnetic fields write: 

0, (V=9 <T <P >lt ) + (v^/2) F«p * v (<p) = 0, %(<p) = dV/dtp, (2.2) 

d u (V=g^y&)) = o. (2.3) 

In accordance with the assumption, made above, the equation (2.3) is easily inte- 
grated at r > 0: 



;F 01 (r) = qP(ip)/J-g' = qP((p)y/l - kr 2 /R 3 r 2 , (2.4) 



R 6 r 4 



where q = const, P(tp) = l/tf(<p) and -g' = -g/sm 9 -- (1 _ fer2) - 

The equation (2.2) for </?(r) in this case metamorphoses to the equation with "in- 
duced nonlinearity" [5]: 

V /Z ?( \Fd'9 11 <P')' = Q 2 9u P v . (2.5) 

Let us make the following assuption of the cosmological character of time in Robertson- 
Walker Universe. Let us suppose that the cosmological time scale is much greater 
than the usual time scale. In other words in the case considered R(t) can be interpreted 
as a constant. Then it is also easy to find the first integral and solution in quadrature 
for the equation (2.5): 

<p' = -V2qVP +C/Rr 2 y/1 - kr 2 , C = const, (2.6) 

J d<p/VP + C = Vl - kr 2 /Rr + C 3 . (2.7) 

The regularity condition of Tq at the center leads to the fact that C — 0. Choosing 
P(<p) in the form 

P{tp) = 1/V(<p) = J 2 - A ' a (l - J 2 '^ , (2.8) 
where J — Xip, a — 2n + 1, n — 1, 2 • • •, for ip(r) one gets: 



(f(r) = - 



1 — exp 



Ra 

It is obvious that at r — > ip(0) — > 1/A, and beginning with some 



(2.9) 



r = r c = 2V2q\/^J(R 2 a 2 Ci + 8kq 2 \ 2 ), 

ip(r) becomes totally imaginary as in this case the square bracket possesses negative 
value. As far as we are dealing with a real scalar field, ip(r) at r > r c becomes non- 
physical. So without losing the generality we may write that at r — > r c , <p(r c ) — > 0. 
(An illustration of the inverse interaction function P(v) and the scalar field obtained 
is given in Figure 1 and Figure 2.) 
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Let us write the energy-momentum tensor for the interacting fields: 

t; = - F^^iy) - s;c. (2.10) 

From (2.10) we find the density of field energy of the system: 

_ 3 q 2 P 

and total energy 

Thus, we came to the conclusion that energy density Tq and total energy of the configu- 
rations obtained do not depend on the conventional values of the parameter k = 0, ±1. 

3 Stability problem 

To study the stability of the configurations obtained we will write the linearized equa- 
tions for the radial perturbations of scalar field. Assuming that 

<p(r,t) = <p(r) + £(r,t), £ < ^, (3.1) 

from (2.2) in view of (2.5) we get the equation for £(r, t) : 

P + 3 — P - 1 - P" - 2 - 3kr2 P' + q2Pw P - C3 2) 

4 + 6 R^ R 2 4 tR 2 « + R 4 r 4 « - U - V-V 

As far as according to the assumption the external gravitational field is cosmological 
one, we can consider that R(t) is a slowly varying time-function: R(t) ~ 0. Assuming 
that 

£(r,t) & v(r)exp(-iQt), Q = uj/R, (3.3) 

from (3.2) we obtain 

(1 - kr 2 )v" + (2/r - 3kr)v' + {uo 2 - ^ E t)v = 0. (3.4) 

Let us first consider the case when k = +1. Then substituting v(r) = y(x), where 
x = 1 — 1/r 2 , from (3.4) we get the equation 

2 2 7~~> 

4*^ + 2y s + (—^— 2 - q -^F) V = °> ( 3 - 5 ) 



;i - x) 2 r 4 

which for y(x) = u(z), x = z 2 , takes the form 



Further substitution 

77(C) = u(z)/VT^, 2 = -the, 
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leads the equation (3.6) to the normal form of Liouville [11] 

ifc C + {u 2 -I- t^e. S ech 4 () V = 0. (3.7) 

In case of A; = —1 the equation (3.4) can analogously be transferred to the form (3.7) 
doing the following substitutions: x = 1 + 1/r 2 , x = z 2 and z = th(. 
At last in case of k — from (3.4) we get 

W" + (u 2 - ^) W = 0, (3.8) 

where W(r) = r • v(r). 

Using the form of P w from (2.8), we come to the conclusion that for a > 5 the 
expressions of the potentials 



\4M = 1 + ^sech 4 C and V^) = Rl 



r • 



tends to +oo at r — > as well as at r — > r r = 2y^<? a = ^ ^ means ^hat 

for a > 5 for P(<p) given by (2.8) the configuration obtained is stable for the class of 
perturbation, vanishing at r = and r = r c . 

In stability can be assured in general introducing the variable 



c = - 



V 

dr If dr 



V^g 11 RJ rVl - kr 2 
and rewriting the equation for perturbation in the form 

+ (fi 2 -g 2 P w K = 0. 

The equation mentioned possesses at Q = nonnegative solution £ = — dip/ d(, which 
according to the Sturm theorem corresponds to the absence of "coupled" state with 
tt 2 < 0. 



4 Conclusion 

Thus, we obtain the object with sharp boundary, described by the regular function 
tp(r). In the center of the system r = y?(0) — > 1/A, and at some critical value 
of radius r — r c function ip(r) possesses trivial value. The configuration obtained, 
possesses limited energy density and finite total energy. From (2.12) it is explicit that 
the expression for energy does not contain r, defining the size of droplet. It means that 
the droplets of different linear sizes up to the soliton with r c — > oo share one and the 
same total energy. For different values of k, the field function ip(r) changes it's form. It 
is noteworthy to notice that at r c — > oo for k = droplet transfers to usual solitonian 
solution, while in case of k — ±1 this type of transition remains absent. It should 
also be emphasized that the values k — ±1 enforce the stability of the configurations 
obtained, which is obvious from the expressions of V (ip) and V± (</?). 
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